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Abstract 

A field-theoretic parametrization is proposed for nucleon electromagnetic form factors at mo- 
mentum transfer less than 600 MeV. The parametrization is part of a larger effective field theory 
lagrangian that is Lorentz covariant and chiral symmetric, and that has been used to successfully 
describe bulk and single-particle properties of medium to heavy mass nuclei. The parametrization 
is based on vector meson dominance and a derivative expansion of nucleon couplings to the elec- 
tromagnetic fields. At lowest order in the expansion, it is possible to fit all four parameters to 
modern data on the rms radii of the nucleon form factors. At next-to-leading order it is possible 
to fit the form factors to within a few percent up to momentum transfers of 600 MeV. The vec- 
tor meson dominance contributions are crucial in this fit, since a simple expansion in powers of 
momentum transfer would require many, many terms to achieve comparable accuracy. The ability 
to fit single-nucleon form factors up to 600 MeV momentum transfer makes possible the study of 
two-body electromagnetic exchange currents within this effective field theory framework. 
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I. INTRODUCTION 



Lorentz-covariant meson-baryon effective field theories of the nuclear many-body prob- 
lem (often called quantum hadrodynamics or QHD) have been known for many years to 
provide a realistic description of the bulk properties of nuclear matter and heavy nuclei. 



(For reviews, see Refs. Ill 
been proposed 0, H, lol. llO. 
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5|,16|.) Recently, a QHD effective field theory (EFT) has 
12l | that includes all the relevant symmetries of the under- 
lying QCD. In particular, the spontaneously broken SU{2)l x SU{2)ji chiral symmetry is 
realized nonlinear ly. The motivation for this EFT and some calculated results are discussed 

in Refs. 0, 0, Q M, M, M, H H E, H in . 

This QHD EFT has three desirable features: (1) It uses the same degrees of freedom 
to describe the currents and the strong-interaction dynamics; (2) It respects the same in- 
ternal symmetries, both discrete and continuous, as the underlying QCD (before and after 
electromagnetic interactions are included); and (3) Its parameters can be calibrated using 
strong-interaction phenomena, like vrN scattering and the empirical properties of finite nuclei 
(as opposed to electroweak interactions with nuclei). It thus provides a natural framework, 
based on a single lagrangian, for discussing the roles of one-body and two-body currents in 
nuclear electromagnetic interactions. 

The nucleon electromagnetic (EM) structure (form factors)js described in this EFT using 
a combination of vector meson dominance (VMD) 0, and a derivative 

expansion for nucleon interactions with the EM field. In the applications of this EFT to 
nuclear structure noted above, however, only the lowest-order derivative couplings were 
included, so that the form factors provided an accurate description of the single-nucleon 
electron scattering data only up to roughly 250 MeV momentum transfer. In contrast, if 
one is to study two-body (exchange) currents, one must reproduce the single-nucleon form 
factors accurately up to at least 600 MeV momentum transfer, where two-body contributions 
are expected to be visible. This momentum scale should be accessible in this low-energy 
hadronic EFT 0, H, Hsj 



Our motivation for this study of nucleon form factors is twofold. First, we want to update 
the lowest-order fits of Ref. p] to include the large amount of low-energy, high-precision data 
that became available in the early 2000's. Second, we extend the fits to the next order in 
momentum transfer and show that the form factors will accurately reproduce the empirical 
results up to roughly 600 MeV momentum transfer. This will make them suitable for studies 
of exchange currents within the QHD EFT. 

In the past ten or fifteen years, much new data on the nucleon EM form factors have been 
obtained using both unpolarized electron scattering and polarization transfer. (For recent 
reviews, see Refs. j29|, [sO].) There have also been numerous attempts at fitting the improved 



data set; for example, see Refs. [3l|, |32|, |33j. For the present study, we are most interested in 
the work of Kelly who achieved excellent fits with a small number of free parameters. 
In particular, the fits are good enough over the momentum transfer range of interest to us 
that we will simply fit our EFT parameters to Kelly's analytic results rather than to the 
data itself. Since our best fits reproduce Kelly's at the few percent level, this procedure is 
justified. 

One of our interesting results is that a straightforward expansion of Kelly's analytic 
results is inadequate unless many, many terms are retained. (Here = —q^ is the square 
of the spacelike four- momentum transfer.) The presence of the VMD contributions in the 
EFT approach greatly improves the situation. Moreover, it is important to include the new 
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EFT parameters in such a way that the error is minimized for the whole relevant range of 
not just ^ 0. 

This paper is organized as follows: In Sec. 11, we return to the lowest-order parametriza- 
tion of Ref. 0] and re- fit the EFT parameters to the new data set. This allows us to 
determine mean-square radii for all four form factors (neutron/proton-electric/magnetic). 
In Sec. Ill, we extend the EFT lagrangian by introducing new parameters and use them to 
fit the higher-momentum transfer behavior of the form factors. Sec. IV is a brief Summary. 



II. RE-FIT OF LOWEST-ORDER PARAMETERS 

In this section, we consider the form factors as described in Ref. 0]. We follow the 
conventions of Refs. 0, [l2|. Rather than work with the Dirac (Fi) and Pauli {F2) form 
factors, defined in terms of the nucleon EM vertex as 

r'^ = Fi(QV + i^2(Q^)^, (1) 

where M is the nucleon mass and F2 contains the anomalous magnetic moment, here we 
will primarily concentrate on the Sachs form factors 

= - TF^m , = F,m + F^m , (2) 

where r = Q^/4M^ = — g^/4M^ in terms of the four-momentum g'^, and we have not 
distinguished the charge states. The charge states are written in terms of the isoscalar (0) 
and isovector (1) parts as, for example, 

Fp = + , F„ = - . (3) 

The simple parametrizations used by Kelly j32| take the form 



En h 

1 + Tr=\ hr^ 



Gm = . , (4) 



which guarantees the correct asymptotic dependence at large Q^: G{Q^) oc Q~^. This will 
not concern us, as we are interested in parametrizing the form factors at small Q^. With 
n = 1 and Oq = 1, this parametrization gives excellent fits to Gep, GMn/jip, and GmuI l^n 
(where /ij is the full magnetic moment) using four parameters each |32|. For Geui Kelly 
follows the so-called Galster parametrization 34l |: 

At 

GEnm = Y:^Gom . (5) 

where the dipole form factor is 

and A and B are fitted parameters. 

For our parametrization, we use set Q2 of Ref. 0] . This provides an accurate fit to bulk 
and single-particle nuclear properties and leaves the anomalous coupling to the isoscalar 
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TABLE I: Coupling parameters from set Q2 [7|]. Note that the nucleon couphngs to the omega 
and rho mesons {g^ and gp) are determined from the empirical properties of nuclei. 



/5(o) 




9v 


9p fv 


fp 


0.01181 


-0.1847 


12.2148 


8.5572 


4.264 



vector meson ("omega") undetermined; we will determine it here. We will need the mass 
parameters 



M 


= 939 MeV 


= 4.7585 fm"^ 


my 


= 782 MeV 


= 3.963 fm"^ , 


nip 


= 770 MeV 


= 3.902 fm"^ , 



(7) 



the anomalous magnetic moments 

A„ = 1.793 , A, 



-1.913 



(8) 



the couplings in Table I, and the electromagnetic coupling g.y = 5.0133, which follows from 
the decay width TpO^e+e- = 6.8 keV. 

We now proceed with the fits to the new data. The calculations use the lagrangian of 
Ref. 0] and are performed at tree level. Starting with the proton electric form factor, 



Gep{Q^ 



.(0) , fj,.(0) 



l_(;5(0)+^(i) 



9p_ 



2M2 



, fpQ'\ Q' , Q' 



4M2 J Q^ + ml 



If we define the mean-square radius as 



^Ep = -6 



dGEpiQ' 



dg2 



Q2 = 



(9) 



(10) 



then 



Ep 



+ 



M2 Sg^ml 



+ 6 



+ 



9p 



M2 g-ymj, 



3Ap 
2M2 



3Ap 
2M2 



(11) 



where the mean-square radii on the right-hand side are the isoscalar and isovector values for 
the Dirac form factor Fi. Inserting the Q2 parameters leads to 



Ep) 



\l/2 



0.862 fm 



(12) 
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which also agrees with the result in Ref. 13 • 
Turning now to Geu, we have 



.(0) 



4M2 



2M2 3g^ 



4:M^ J Q^ + ml 



so that 



^2 
'En 



1 
2 



M2 3^^m2 



-6 
3A 



+ mj 



An 



4M2 



(13) 



+ 



3A„ 
2M2 



(14) 



If we set (r2)(o)i = (r^)^]^)]^, as in Ref. [7|, we then find r|;„ = — 0.127fm^, in significant 
disagreement with Kelly's value of —0.112 ± 0.003 fm^. We conclude that the new data 
shows that 

(^'){0)i-(O(i)i = 0.0294fm2. (15) 

With this information, together with Eq. ffTTl) . we can determine two distinct radii for the 
Dirac form factor: 



(r2) = 0.799 fm , (r^) J/)\ = 0.780 fm 

in contrast to the assumption made in Ref. Q. 

We now consider the magnetic form factors, beginning with 



(16) 



l-(/3(°)+/5«: 



2M2 
9pil + fp) Q' 



2g, Q^ + mj 



gyjl + fv) Q' 



(17) 



As expected, for — > 0, GMpiQ"^) — > 1 + Ap = yUp. Thus we normalize Gmp by dividing by 
/ip, leading to the mean-square radius 



1 



Mp 



6 



2(1 + Ap 
^ 2(T+A, 



M2 ?>g^ml 



6 



M2 (7^777,2 



p/ J 



2 1 + A, 



g^m^^ g^ml J 
((r')(o)i + + (Ap + Xn){r^)m + (Ap - Xn){r^)(i)2) 



0.7191 fm^ 



(18) 
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where the numerical value is taken from Kelly. Note that the magnetic radii depend on both 
the Dirac and Pauli mean-square radii. 
For the neutron, 



GMniQ ) — -^1 



(0) 



+ 



2M2 
9p{l + fp) Q' 



'29, 



The mean square radius is 
1 



Mn 



1\r. 



6 



+ 



M2 ^ 



M2 



3gy + ml 



9p 



(19) 



p/ J 



2Xn yg^ym. 



Qfp9p \ 
9,rnl) 



2Xr 



0.8226 W 



(20) 



With the results in Eqs. ffTB]) . ffTSl) . and fl2U]) . it is now simple algebra to determine the 
rms radii for the Pauli form factor, with the results 



/ 2\l/2 



1.30 fm 



/ 2\l/2 



0.896 fm . 



'(0)2 ~ ^^^^ ' y /(i)2 — ^•"■^^ ^^^^ • (21) 

The isovector radius is consistent with Refs. 0, [35| . At this time, we have no way to estimate 
errors for any of these derived radii. 

Equation (1201) also shows the relationships between the Dirac and Pauli rms radii and the 
C((5^) parameters in the EM coupling expansion (i.e., /J'-^-', /„, /p). Simple inversion of 
these results allows us to present parameters determined from the new single-nucleon data 
with the nuclear couplings of set Q2: 

/5(o) = 0.0670 , = -0.2402 , = -0.3279 , fp = 4.4198 . (22) 



III. INCLUSION OF HIGHER-ORDER COUPLINGS 

So far, we have adjusted the parameters that enter at 0{t) to achieve the modern values 
for the four rms radii. Now we want to go to higher order in r and see what adjustments are 
necessary to reproduce the Sachs form factors up to Q ~ 600 MeV or r ~ 0.1. We can relate 
these adjustments to higher-order terms in the derivative expansion of the EM lagrangian. 

For the purposes of this work, we take the nucleon part of the EM lagrangian as 

£em = -eApN^lil + rs)N - -^F,,N\a^''N - ^{d'^ F,,)NprN 

- ^{d,d^'F,,)NX'a>^'^N - ^^{d'd''F,,)N(3YN 

- 7^,id'd.d'F,,)N\"a>^^N + ■ ■ • , (23) 
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where A^^ and F^'^ are the electromagnetic fields. All of the constants A, /3, A', A" have 
the isospin structure 

1 111 

A = Ap-(1 + rg) + A„-(l - rs) = -{Xp + AO + ^r^iXp - A„) 

= A(°) + A^^Va , etc. (24) 

As shown in Ref. [13], the isovector parts of these constants should be modified to include 
pion interactions to maintain the residual chiral symmetry in the lagrangian. When applied 
to the nucleon form factors, however, these pions appear only in loops, which we are not 
considering, so we have omitted these terms. 

We have already included the A and P parameters in the preceding section. It is easy to 
see that the A' constants enter the magnetic rms radii at the same order in as the vector 
meson couplings and fp. Thus the A' parameters are redundant in our approach and will 
not be considered in the sequel. Thus the four new adjustable constants at our disposal are 
contained in and A". 

It is straightforward to work out the Feynman rules for the new vertices and to construct 
the tree-level contributions to the form factors. One finds that these constants enter the 
Dirac and Pauli form factors at 0{Q'^). Thus they will enter the Sachs form factors at 0{t'^) 
and 0(r3). 

Our strategy is the following: We begin with the magnetic form factors, where both new 
terms enter at C(r^). We numerically adjust the coefficient of this term until we get a good 
fit to the form factor up to r = 0.1. Here we define a good fit as one that minimizes the 
maximum deviation between the fit and the data (i.e., Kelly's results 32|) throughout the 
whole interval < r < 0.1. The result provides two constraints between the (3' and A". We 
then turn to the electric form factors, where the new terms enter at both O(t^) and 0{t^). 
The coefficients of both of these terms are then adjusted, consistent with the constraint, 
until a good fit is obtained for < r < 0.1. This procedure provides us with four numbers 
that can be used to determine the (3' and A". 

Denoting the new contributions to the Sachs form factors as 5G, we find 

5GMp{r) = 16(A"(°) + A"(i) - - V , (25) 

5GMn{r) = 16(A"(°) - A"(^) - + p'^'^y . (26) 

Taking the best fit to the data leads to Figs. [T] and [21 and the constraints 

a;'-/?; = 2.255, a:-/?; = -2.786. (27) 

The fits are accurate to a few percent. 

For the electric form factors, the new terms are 

SGEpiT) = -16(/?'(°) + - 64(A"(°) + A"(i V , (28) 

6GEn{r) = -16(/3'(°) - (3'^''>y - 64(A"(°) - A"^^ V (29) 
Taking the best fit subject to the constraints above produces Figs. [3] and HI and the values 

p'^ = -0.9438 , = 0.4688 . (30) 

Combining these results with the constraints determines all four new parameters. 

Although the relative error in Geu is as large as 30%, Geh is small in the momentum- 
transfer range of interest. A more relevant comparison is given in Fig. [5l where all four of 
the empirical form factors are shown along with fits to Geu and Gmu- 
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FIG. 1: Proton magnetic form factor as a function of r. The curves represent the data (sohd), the 
lowest-order fit (dashed), and the new fit (dot-dashed). 



IV. SUMMARY 

In this paper we studied a field-theoretic parametrization of single-nucleon EM form 
factors at low energy. This parametrization is part of a Lorentz-covariant, chiral invariant, 
hadronic effective field theory that was proposed to study the nuclear many-body problem. 
We thus have a single lagrangian that describes the nuclear structure, nuclear currents, and 
interaction vertices at low energies. This is a natural framework for discussing the roles of 
one-body and two-body currents in nuclear electromagnetic interactions. 

The parametrization of the form factors is based on a combination of vector meson dom- 
inance and a derivative expansion for nucleon interactions with the EM field. At leading 
order in derivatives, the form factors accurately reproduce the singlc-nuclcon electron scat- 
tering data only up to roughly 250 MeV momentum transfer. To study two-body exchange 
currents, however, one must reproduce the form factors accurately up to at least 600 MeV 
momentum transfer. This paper is proof of principle that by including the next-to-leading 
order (nonredundant) derivatives, one can adequately describe the form factors up to 600 
MeV with our form of parametrization. 
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FIG. 2: Neutron magnetic form factor as a function of r. The curves are identified as in Fig. [H 



We also revisited the leading-order parametrization and used a modern data set to eval- 
uate the expansion coefficients. We found that it is now possible to determine all four rms 
radii for the Dirac and Pauli, isoscalar and isovector form factors, and thus determine all 
four leading-order coefficients from the single-nucleon data, unlike in Ref. 

It is interesting that a simple power-series expansion in powers of the momentum transfer 
squared would require many, many terms to reproduce the form factors up to the desired 600 
MeV. (This is easy to see using ten minutes of simple calculations on mathematica with 
Kelly's [ill fits to the data.) Thus the vector meson dominance contributions are critical in 
allowing us to parametrize the desired data using only the next-to-leading order derivative 
terms. 

While there is much work still to be done before two-body currents can be studied numer- 
ically within this effective field theory framework, we have shown that the one-body current 
can be adequately parametrized to make a study of two-body currents possible. 
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FIG. 3: Proton electric form factor as a function of r. The curves are identified as in Fig. [H 
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